Regular and chaotic Bose-Einstein condensate in an accelerated Wannier-Stark lattice 
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We investigate a Bose-Einstein condensate held in a quasi-one-dimensional Wannier-Stark lattice 
which is a combination of linear potential with an accelerated optical lattice. It is demonstrated 
that the system can be reduced to a periodically driven Gross-Pitaevskii one, in which we find 
the first exact analytical solution and the regular and chaotic numerical solutions with accelerated 
atomic flow densities. The results suggest an experimental scheme for generating and controlling 
the accelerating regular and chaotic matter-waves. 
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I. INTRODUCTION 



Recently, experimental investigation on the Bose- 
Einstein condensates (BECs) in a Wannier-Stark (WS) 
system have attracted much attention In some ex- 
periments, a BEG of rubidium atoms was created, which 
then was loaded into the optical lattice potential of a 
standing laser light and was adopted to observe the 
collective tunneling effects The optical potentials 

have been used in many theoretical and experimental 
studies of quantum dynamics, for example, the coherent 
pulse output from BEGs in WS system [sj, the observa- 
tion of Bloch oscillations both with single atoms and with 
a BEG in an accelerated standing wave 0, [J , the stud- 
ies of WS ladders for the accelerated optical potential 
0. An atomic-scale analog of the kicked rotor is real- 
ized by placing laser-cooled atoms in a pulsed standing 
wave Q. Very recently, M. Gliick et al. investigated the 
properties of a coherent superposition of WS resonances 
Q and gave the lifetime of WS states (3- Ghaotic be- 
haviors have also been found in the WS systems without 
acceleration H, [1, [ld| • It is a natural motivation for us 
to demonstrate the chaotic and regular features for the 
BEG in an accelerated WS lattice. 

On the other hand, it is well-known that the BEG 
governed by a Gross-Pitaevskii equation (GPE) without 
external potential is an integrable system and the inte- 
grability could be easily broken by external potentials 
of different forms So previously, only few analyti- 

cal works concern exact solutions of the system, where 
one-dimensional (ID) stationary systems with some sim- 
ple potentials are treated, such as the infinite or finite 
square- wells [H, [H, [l^, [isl. the step-potentials ,16], 5 
or 5 comb potentials |l7l. [ig. [l9t. linear ramp potential 
[20l . [2l| and optical lattice potentials [H, HI]. Under 
some rigorous conditions on the interaction intensities 
or external potentials, several exact nonstationary-state 
solutions were constructed |2l, j2&| . including the exact 
soliton solutions [1^, [13, H^Tl^rsO]- It is worth not- 
ing that the balance between nonlinearity and dispersion 
was found in the seminal work of soliton |3l| , and the new 
balances between the atom-atom interaction and the ex- 
ternal potentials are demonstrated recently [1^ [13, [l^ . 



By using the balance conditions, although some exact so- 
lutions have been constructed for the GPE with periodic 
potential, however, any exact solution in the noninte- 
grable WS system with a combination of the linear and 
periodic potentials has not been reported yet. 

The aim of this paper is to present the first exact an- 
alytical solution with balance condition and to illustrate 
the regular and chaotic numerical solutions of the accel- 
erated WS system. The corresponding atomic ffow densi- 
ties accelerated by the constant force are demonstrated. 
Based on the relations between the system parameters 
and the solution behaviors, we suggest an experimental 
method for controlling the regular and chaotic states by 
applying the accelerated optical potential and adjusting 
the system parameters. 



II. SIMPLIFICATION OF THE 
WANNIER-STARK SYSTEM 

The mean-field theory is a successful one for describ- 
ing the BEGs. In this theory, the dynamical behaviors 
are governed by the GPE [sj, [1^ , which provides us a 
nonlinear macroscopic quantum system. Let us consider 
a BEG trapped in one-dimensional tilted optical lattice 
potential 



V{x', t') = cos(2/cL^') 4- Fx', 



S.' = x' + \at'^- 



(1) 



here x' and t' are spatial and time coordinates, 
Vb cos(2fcL^') is the accelerating optical potential [1, [1| 
with strength Vq, wave vector fci and acceleration a, and 
F = ma is a constant force 0] with m being the atomic 
mass. Due to this force, a "tilted" potential is produced 
that leads the atoms to accelerate in x direction with lin- 
early increasing flow density and makes the atoms tunnel 
out of the optical traps. The corresponding dimension- 
less GPE reads as 



9-0 



dx^ 



+ [V„cos{20 + ax + gu\^\'']^, (2) 
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where the dimensionless spatial and time coordinates 
are x = kj^x' and t — Erf /h. The wave function ip 
has been normalized in units of the radical VAJl, the 
potential depth Vo is normalized by the recoil energy 
Er = li^k\/{2m). The constant force is rescaled by the 
unit k^Er, so all variables and parameters in equation 
(3) are dimensionless. In such units, the interatomic in- 
teraction intensity related to the s-wave scattering length 
as is in the form gid = Aas/{kLl'^) with Ir = \J h/ (mcor) 
being the radial length of harmonic oscillator. 

In order to get a simple description and better under- 
standing of BEC dynamics, we let the wave function be 
in the form 



(3) 



where undetermined function t) may be real or com- 
plex, which is normalized to the total number of atoms, 
/ v?{£^,t)dx = N. Given Eq. (2), we perform the calcu- 
lations 



9V 
dx'^ 



\dt 



+ 2at^')e-'("-*+^°'*') + {ax + a't')i,, 
at, J 



2at 

9^ " V9^"*^"*9^y 



du\ 

e 



ah^ip. (4) 



Substituting Eq. (4) into Eq. (2) yields 



.du d^u . ^ , 

= - ^ + [^^ocos(20 + gid\u\]u 



(5) 



in which the linear potential is removed and its effect is 
included in the parameter a. When the optical potential 
is switched off, Vq = 0, Eq. (5) becomes a standard non- 
linear Schrodinger equation (NLSE), whose single-soliton 
and multisoliton solutions are well-known for us. It is 
worth noting that in the transformation x x + wk for 
n = 0, 1, 2, • • • , Eq. (5) and its solution t) are kept 
such that Eq. (3) gives n solutions 



IpniXj t) = ijj{x + riTT, t) = t)e 



(6) 



The solutions with different n possess the different phases 
and the same amplitude. 



III. EXACT PERIODIC WAVE WITH 
ACCELERATED FLOW 



where /i is a constant adjusted by the normalization con- 
dition and can be call the chemical potential. By using 
the balance technique, we establish the balance condition 



9id\ 



(9) 



and then reduce Eq. (7) to the linear Schrodinger equa- 
tion 



(10) 



The exact solution of Eq. (8) must obey Eqs. (9) and 
(10) simultaneously. General solution of the complex 
equation (10) can be written as 



= [A + iC] cos^ + {B + iD)ii\i\^, 



(11) 



where A, B, C and D are real constants which are de- 
termined partly by the balance condition (9). From Eq. 
(11) we construct the quadratic norm 



101^ 



A^ + + {B^ + 
(AB + CD) sin{2^). 



A^~C^)sm^^ 



(12) 



Comparing Eq. (9) with Eq. (12) and noticing cos(2^) — 
1 — 2 sin^ ^ produce the algebraical equations 

gid{A^ + C^)=li-Vo-l, AB + CD = 0, 
gid{B^ + D^-A^-C^)^Vo, (13) 

which denote a group of indefinite equations with infinite 
numbers of solutions. 

The existence of multiple solutions of Eq. (13) implies 
that phase of the exact solution 



arctan 



{Acoa£_ + Bsin^) ' 



(14) 



has some arbitrariness. However, quadratic norm of 
the solution (7) is determined uniquely by Eq. (9), 
= [/i — 1 — Vq cos{2£^)]/gid- This means that un- 
der the exact state the atomic density profile shapes the 
periodic wave-packets which propagate with acceleration 
a. Therefore, we can achieve the accelerated transport of 
BEC, through the considered exact solution. Application 
of the normalization integral yields the average number 
of condensed atoms per well [22] 



We are interested in the exact analytical solution and 
regular and chaotic numerical solutions of Eq. (5). Notic- 
ing that ^ and t in Eq. (5) are two independent variables, 
we can rewrite the function u(^, t) in the separation form 
of variables 



(7) 



and transform Eq. (5) to the ordinary differential equa- 
tion 



li<j> - 



+ [yocos(20+ffid|</'|']<^, 



(8) 



N' = (tt)- 



\c^{x)\^dx = {ii-l)/g 



Id, 



(15) 



which determines the chemical potential as /i = gidN' +1. 
Given the chemical potential, the exact atomic density 
reads 



= i?2 = Ar'--^cos(2e). 
9id 



(16) 



The result is accurately adjusted by the experimental 
parameters iV', Vq, gid, and wave vector fc^ and accel- 
eration a implied in ^. 
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Clearly, Eq. (8) is a periodically driven GPE of the 
spatiotemporal evolution in which the well-known Smale- 
horseshoe chaos exists for a certain parameter region [s^ . 
Writing the solution (p in the form of (f) = R{(.) exp[i6'(^)] 
and inserting it into Eq. (8) lead to two couple equations 



= R 



gidR"" + [Vocos{2^)~ ^l]R, (17) 
0. (18) 



The square of the modulus \tp\'^ = |</)p = R^ denotes 
atomic number density and the total phase reads 

Q{x, t) = e{S) + ann)t - {axt + -a'^t^). (19) 

3 

The both are associated with the velocity field v and 
flow density J, through the formulas v = hOx/m and 
J = vR^. Integrating Eq. (18) yields the part phase 



So the velocity field and flow density become 

Jo 



(20) 



v{^,t) = (%-a<) = 

m m 

m 



R'iO 



at 



(21) 



where Jo — S^{^o)R'^{£,o) denotes an integration constant 
determined by the flow density at the point S,o ^ xq+ at^ 
for the initial time to and boundary coordinate xq. Sub- 
stituting Eq. (16) into Eq. (21) produces the exact flow 
velocity and flow density. The linear term of time in flow 
velocity implies the BEC superfluid being accelerated. 



IV. REGULAR AND CHAOTIC NUMERICAL 
SOLUTIONS 

In general, the balance condition (9) cannot be satisfied 
such that we have to solve Eq. (17) for the modulus R. 
Applying 6^ = Jo/i?^ of Eq. (20) to Eq. (17), we arrive 
at the decoupled equation 



7^2 

•^0 



[Vo cos(20 - fAR. 



(22) 



This is a real equation on the accelerated reference frame, 
which can be reduced to the parametrically driven Duff- 
ing equation [sl, [11] for the case Jo = 0. The Smale- 
horseshoe chaos in Jo = case has been widely inves- 
tigated. Generally, Jo ^ 0, the system becomes more 
complicated and the chaotic property may be kept. Un- 
der balance condition (9) we can prove directly that Eq. 
(16) is an exact special solution of Eq. (22) for a particu- 
lar integral constant. When the periodic driving is weak 



enough, the chaotic perturbed solution has been con- 
structed for the system without linear potential [H, [3^ . 
For most of parameters and initial data Eq. (22) is not 
analytically solvable, which necessitates the numerical so- 
lutions. 

At the initial and boundary point ^o , the values of the 
velocity field and number density can't be determined ex- 
actly in experiment due to the fluctuation of the atomic 
thermal cloud, we can randomly choose some possible 
values to perform the numerical calculations. In order 
to explore the analytically insolvable system (22), we nu- 
merically solve it by using the MATHEMATICA code 

T — tt; e[{Rnew-,vnew-}] := {R[T], v[T]}/ .Flatten 
[NDSolve[{R'[^] v[^],v'[^] g,,R[^]^ - f,R[^] 
+Jq/R[^? + VoCos[2£,] * R[£.],R[0] == Rnew, 
v[0] =~ vnew}, {R, v}, {^, 0, T}]]; Do[piCi — ListPlot 
[Drop[Nestlist[e, {Random[Real , {—0.5,0.5}], 
Random[Real, {-0.5, 0.5}]}, 5100], 100], {i, 1, 10}] 

to make 10 groups of orbits on the Poincare section of 
the equivalent phase space (i?, R{) for the random initial 
conditions {i?(Co) € [-0.5,0.5], i?^(^o) e [-0.5, 0.5]} and 
different parameter sets. Each of the groups contains 10 
orbits, which corresponds to one of the following cases: 

Case 1. For the parameter set gia = —1, /i = 
-0.5, Jo = 0.01, Va = 0.05 afl the 10 orbits are reg- 
ular, whose 3 typical profiles are shown in Fig. 1. 

Case 2. By increasing the strength of lattice potential 
to Vb ~ 0.2 and keeping the other parameters as in case 
1, we observe 3 regular orbits of Fig. 1 and 7 chaotic 
orbits, whose 2 typical profiles are shown in Fig. 2. 

Case 3. After further increasing the strength of lattice 
potential to Vb = 0.5, all the orbits become chaotic. 

Case 4. By increasing the fiow density to Jo — 0.16 
and keeping the other parameters as in case 3, we find 
that all the orbits become regular. 

Case 5. By increasing the strength to Vb = 3 and 
keeping the other parameters as in case 4, we observe 4 
regular and 6 chaotic orbits. 

Case 6. After further increasing the strength to Vb = 5, 
all the orbits become chaotic. 

Case 7. For the parameter set gid — —1, /i — 0.5, Jo — 
0.16, Vb — 0.5 with positive chemical potential /i all the 
chaotic orbits in case 3 are transformed to regular ones. 

Case 8. By increasing the strength to Vb = 2 and 
keeping the other parameters as in case 7, we observe 2 
regular and 8 chaotic orbits. 

Case 9. After further increasing the strength to Vb = 4 
and keeping the other parameters as in case 8, all the 
orbits become chaotic. 

Case 10. For the different parameter sets with posi- 
tive interaction gid we make many regular orbits and no 
chaotic orbit is found. 

The above results imply that for the attractive inter- 
action and negative chemical potential increasing the lat- 
tice strength can strengthen chaoticity of the system. 
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FIG. 2: Typical chaotic orbits on the Poincare sections of the dimensionless 'phase space' {R, R^) for ten different initial 
conditions and/or parameter sets. Here the phase orbits evolve in a finite region and exhibit the confused state distributions. 



Conversely, increasing the flow density can weaken the 
chaoticity. On the other hand, after changing the chem- 
ical potential from negative to positive, the chaoticity is 
weakened cfl'ectivcly. Finally, for the positive interaction 
gid we have not found the chaotic orbits. The 4 typi- 
cal regular orbits and 4 typical chaotic orbits among the 
10 X 10 = 100 orbits are shown in Fig. 1 and Fig. 2 
respectively. They are associated with different initial 
conditions and/or parameter sets. As the examples of 
densities of atomic number we plot their spatiotemporal 
evolutions for the parameters of Case 2 and two fixed ini- 
tial conditions as in Fig. 3. Figure 3(a) corresponds to 
the first closed orbit of Fig. 1 and describes a quasiperi- 
odic evolution thereby. Figure 3(b) is associated with the 
first chaotic orbit of Fig. 2, which possesses obvious ape- 
riodic! ty. From Eq. (21) the corresponding flow densities 
can be easily illustrated. They will increase linearly in 
time that leads the condensed atoms to tunnel out of the 
optical traps. These results reveal the relations between 
the system paramcitcrs and the solution behaviors, and 
suggest a method for controlling the regular and chaotic 
states. 




20 40 

(b) ^ 



FIG. 3: Spatiotemporal evolutions of the densities of atomic 
number for the parameters of Case 2 and the initial conditions 

(a) i?2(0) = -0.07793183579, i?|(0) = -0.199975080313, (b) 
7?2(o) = -0.48451144892, R^O) = 0.31792019031. 



V. CONCLUSIONS 

In summary, we have investigated a BEC interacting 
with an accelerated WS potential. Using the mean-field 
method and the macroscopic one-body wave function, we 
seek the exact analytical solution and the rc;gular and 
chaotic numerical solutions of the system. It is demon- 
strated that after the linear potential being removed by a 
function transformation, the governing equation becomes 
a periodically driven CPE on the accelerated reference 
frame. With the help of the balance condition, we estab- 
lish the first exact analytical solution of the WS system, 
which is accurately controlled by the experimental pa- 



rameters. Further writing the solution of CPE in the ex- 
ponential form, we obtain the equation of modulus, which 
contains the parametrically driven Duffing equation. The 
well-known Smale-horseshoe chaos and quasiperiodic or- 
bits on the Poincare sections of the dimensionless 'phase 
space' are shown numerically for different initial condi- 
tions and parameter regions of different chaoticity. The 
accelerated atomic flow densities are demonstrated for 
both the regular and chaotic states. 

It is well known that chaos could emerge in the proces- 
sions of BEC collapsing and may play a destructive role 
for the BEC system. Therefore, predicting and control- 
ling chaos are quite important for the creation and appli- 
cation of BEC. Our analytical and numerical results have 
supplied a method for controlling the regular and chaotic 
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states, through the appHcation of accelerated optical po- and 10875039. 
tential and the adjustments of system parameters. 
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